Application of Image Analysis for Characterization a: 

of Spatial Arrangements of Features in Microstructure 7// 
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A number of microstructural processes are sensitive to the spatial arrangements of features in 
microstructure. However, very little attention has been given in the past to the experimental 
measurements of the descriptors of microstructural distance distributions due to the lack of prac- 
tically feasible methods. We present a digital image analysis procedure to estimate the micro- 
structural distance distributions. The application of the technique is demonstrated via estimation 
of K function, radial distribution function, and nearest-neighbor distribution function of hollow 
spherical carbon particulates in a polymer matrix composite, observed in a metallographic section. 
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I. INTRODUCTION 

SPATIAL arrangements of features in a micro- 
structure can be quantified in terms of the distances be- 
tween the corresponding features. Distribution of distances 
between features in a microstructure is an important 
microstructural attribute. The microstructural distance 
distributions affect the microstructural processes that in- 
volve short range, and/or long-range particle-particle 
interactions (or feature to feature interactions, in gen- 
eral). For example, diverse phenomena, such as for- 
mation and growth of microcracks during thermal cycling 
of metal matrix composites, 11 * 21 damage accumulation and 
fracture behavior of composites, 13,41 particle coarsen- 
ing, 15,61 liquid-phase sintering, 18 91 microstructural evo- 
lution during solid-state transformations where nucleation 
occurs on the second-phase particles, 110 1,1 microvoid in- 
duced ductile fracture processes in metals and al- 
loys, 112,131 creep cavitation, 1141 and creep crack growth, 1151 
crucially depend on the distances between the relevant 
microstructural features. 

The microstructural distance distributions may be sig- 
nificantly altered when a material is processed in re- 
duced gravity or microgravity, particularly if a liquid phase 
is involved. In such a case, it is desirable to quantify the 
changes in the microstructural distance distributions due 
to reduced gravity and to quantify and model the role of 
gravity in the microstructural evolution processes. 

Very few investigations deal with the quantitative ef- 
fect of the spatial distribution of microstructural features 
on the physical/mechanical properties of materials or on 
the microstructural evolution during materials process- 
ing. An important reason for this is the lack of well- 
developed practical techniques to quantify the spatial 
distribution of microstructural features and, as a result, 
the lack of realistic and flexible quantitative models to 
represent the spatial distribution in the computer simu- 
lation or analytical theoretical studies. 
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It is the purpose of this article to present an image 
analysis procedure to obtain the basic experimental data 
necessary to quantify the spatial arrangement of micro- 
structural features in metallographic sections. To char- 
acterize spatial distributions, the necessary basic data 
consist of centroid coordinates and sizes of the micro- 
structural features of interest. These raw data can be pro- 
cessed to compute the descriptors of spatial order, such 
as nearest-neighbor distribution, 116 17 181 radial distribu- 
tion, 119 ' 221 K function, 119 ' 221 pair-correlation func- 
tion, 12122 * 231 or parameters such as short range and long 
range microstructural gradients. i24) An image analysis 
procedure was developed to automatically extract cen- 
troid coordinates (and feature size, etc.) of the micro- 
structural features observed on different Fields of view 
and referred to the same origin: this is essential to com- 
pute the distances between the particles that may not be 
in the same Field of view. The image analysis algorithm 
provides a practical procedure to acquire the centroid co- 
ordinates data of thousands of particles to quantify the 
short-range as well as long-range afFinity (or lack of it) 
among the features. A simple computer program utilizes 
these raw data to compute the descriptors of spatial dis- 
tribution, such as nearest-neighbor distribution and ra- 
dial distribution. 

For population of spherical particles, stereological re- 
lationships are available to estimate the three-dimensional 
(3-D) radial distribution from the corresponding experi- 
mentally measured two-dimensional (2-D) radial distri- 
bution of the particle sections observed in a metallographic 
plane. 1,9 21 231 In the case of composites having continu- 
ous aligned fibers, the spatial distance distributions of 
Fiber centers observed in a metallographic plane perpen- 
dicular to the Fibers are sufficient for the characterization 
of the spatial distribution of the fibers and hence no ster- 
eology is necessary in such cases. 

In the present study, the image analysis procedure is 
applied to quantify the spatial distribution of hollow 
spherical carbon particles in a polymer matrix compos- 
ite. The measurements are performed on two specimens 
having different volume fractions of particulates but the 
same size frequency function A brief background on 
microstructural distance distributions is presented in 
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Section II. The image analysis procedure and practical 
example are given in Sections III and IV. 


II. DESCRIPTORS OF MICROSTRUCTURAL 
DISTANCE DISTRIBUTIONS 


Let N v be the average number of particle or feature 
centroids per unit volume of microstructure. Focus on 
one particle whose centroid is at location (X n Y, t Z t ) 
in the microstructural space. Draw a test sphere of 
radius R around this particle having the center at 
(X h Y n Z,)). Let K*(R) be the total number of other par- 
ticle centroids contained in the test sphere. Obviously, 
the quantity JC(R) depends on the location of the test 
sphere (i.e., (X n Y n Z,)), the distance /?, and the spatial 
arrangement of the particle centroids in the micro- 
structure. An average value, K(R)> can be obtained by 
averaging fC(R) over all the particle centroid locations 
(X n Y n Z t ) in the microstructural space. The function K(R) 
is called the K function, ll9 ~ :31 and it gives the average 
number of particles (strictly speaking, particle centroids) 
in a sphere of radius R drawn around the centroid of an 
arbitrary particle or feature of interest in the micro- 
structure. For a collection of randomly distributed par- 
ticles, K(R) is equal to (4 tt7? 3 /3) • N v ; a value higher than 
this signifies clustering, and a lower one indicates 
repulsion. 

The radial distribution function G(R) is an alternate 
descriptor of the spatial arrangement of particle (or fea- 
ture) centroids. It is defined as follows: 119 20 211 


C(/?) = 


1 

(4 77 R 2 N v ) 


dm) ] 

dR J 


[ 1 ] 


The quantity (4ttR 2 * N t dR) • G(R) is equal to the av- 
erage number of particle centroids in a spherical shell of 
radii R and (R + dR) around an arbitrary particle. For a 
random spatial distribution of the panicles, G(R) does 
not depend on R and it is equal to 1. For a clustered 
distribution of panicle centroids, G(R) is greater than 1 
for small values of R. It is composed of a combination 
of delta functions for an ordered array of panicles. 

Note that the radial distribution function is related to 
the derivative of the K function (Eq. [1]), and hence it 
is more sensitive to changes in the spatial order than is 
the K function. The radial distribution function is par- 
ticularly useful to quantify short-range particle-panicle 
interactions. For most of the spatial arrangements of par- 
ticles, G(R) approaches 1 for a sufficiently large value 
of /?, and it approaches zero as R tends to 0. 

The radial distribution function contains information 
concerning the spatial arrangement of particle centers, 
but it gives no information concerning correlations be- 
tween sizes and relative locations of the particles. For 
example, it does not tell us if the small particles are pref- 
erentially located near the larger particles, and so on. 
The pair correlation function contains such information. 
It is defined as follows: ,2, • 22 ' 23, 


or,, r 2 , r) = \ — 2 U ^r,,r,.«) | 

[47r/? 2 • « v (r 2 )J 1 dVydR j 1 J 

n v (f 2 ) is the particle size distribution function, such that 


n v (T 2 ) dT 2 is the number of particles in the size range 
r 2 to (T 2 + dT 2 ) per unit volume of microstructure. The 
value of Q( r,, T 2 , R) is the average number of particle 
centroids belonging only to the particles of size smaller 
T 2 in a sphere of radius R centered on a particle of size 
T,. The value of C(r,, T 2 , R) is the pair correlation func- 
tion for the particles of the sizes T, and T 2 separated by 
distance R . The pair correlation function is 0 for R < 
(r, + T 2 ), and it approaches 1 as R tends to infinity. For 
a random spatial distribution of particle centroids and 
uncorrelated particle sizes, the value of the pair corre- 
lation function is equal to 1 for all the values of /?, T,, 
and T 2 . The pair correlation function is particularly use- 
ful in the study of processes, such as particle coarsening, 
where small particles may be preferentially located near 
the larger ones. 15 61 

In the deformation and fracture processes, such as 
ductile fracture and 112 131 creep crack growth and frac- 
ture, 1151 the nearest-neighbor distribution function and 
average nearest-neighbor distance are the relevant 
microstructural parameters. The nearest-neighbor distri- 
bution function 1 l6 l7J81 is described by the probability 
density function P{R ), such that P(R) • dR is equal to the 
probability that there is no other particle centroid in a 
sphere of radius R around a given particle and there is 
at least one particle centroid in the sphencal shell of radii 
/? and ( R + dR). The average nearest-neighbor distance 
R is given by the following equation: 


R = 



R • P (R) dR 


[3] 


For randomly distributed point particles, R is equal to 
0.554 ' (N v ) 333 . In a manner similar to the nearest- 
neighbor distribution function, one can define 2nd, 3rd, 
. . . nth nearest-neighbor distributions. ,::1 

The materials microstructures are usually character- 
ized from observations on representative 2-D metal Jo- 
graphic sections through the 3-D microstructural space 
of interest. Thus, the basic microstructural data often 
pertain to the 2-D metallographic sections. All the pre- 
ceding descriptors of distance distributions can be anal- 
ogously defined for the 2-D microstructural sections. The 
K function pertaining to 2-D microstructural section, k(r) 
is defined as the average number of particle section cen- 
troids in a circle of radius r, centered on an arbitrary 
particle section. The radial distribution function g(r) per- 
taining to a 2-D microstructural section is defined as 
follows: 


g(r) = 



tt ' r ' N a 



14] 


where N A is the average number of particle section cen- 
troids per unit area of the metallographic section. Note 
that the lower-case symbols denote the distance distri- 
butions in the 2-D sections. The nearest-neighbor dis- 
tribution function in a 2-D section is described by the 
probability density function p(r) % such that p(r) • dr is the 
probability that there is no other particle section centroid 
in a circle of radius r centered on an arbitrary particle 
section of interest and there is at least one particle sec- 
tion centroid in the circular ring of radii r and (r + dr). 
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The pair correlation function for 2-D microstructural 
section c(r,, t 2 , r) can be defined as follows: 


c(r lf x 2 , r ) = 


2'TT'n 0 (r 2 ) 


dq(r u t 2 , r) 


n M (r 2 ) is the particle-section size distribution observed in 
the metallographic section. The value of q(r {y t 2 , r) is 
the average number of particle-section centers belonging 
only to the particle sections smaller than r 2 in a circle 
of radius r centered on an arbitrary particle section of 
size 7 ,. 

If microstructural distance distributions are experi- 
mentally measured on representative 2-D metallographic 
sections, then under certain conditions, the correspond- 
ing 3-D distance distributions can be estimated by using 
stereological relationships. |l9,2 ' 231 For example, if the 
particle shape can be assumed to be spherical, then the 
3-D radial distribution function G(R) and the corre- 
sponding 2-D function g(r) measured on a metallo- 
graphic section are related through the following 
stereological equation: 1191 


8(r ) = {^y} ' J o + m2 )' /2 1 ’ d “ 16] 

where, 



where T* is the average particle radius, and / is the cu- 
mulative size frequency function of the spheres. Note 
that in Eq. [6], the desired 3-D radial distribution func- 
tion G is under the integral sign, and hence numerical 
solution to this integral equation is necessary for the es- 
timation of G. 


III. IMAGE ANALYSIS TECHNIQUE 

To estimate any of the distance distributions in a 2-D 
metallographic section, it is necessary to measure the 
( X , Y) coordinates of the observed particle-section cen- 
ters and the corresponding section sizes. To obtain sta- 
tistically meaningful measures, it is essential to perform 
these measurements on a large number of particle sec- 
tions distributed over many fields of view (or frames), 
and hence automatic digital image analysis is essential. 
for all practical purposes, these measurements cannot be 
performed manually. In the digital image analysis, the 
microstructural distances are measured in units of pixels; 
the absolute sizes are calculated from this information 
by using the calibration factor, which depends on the 
microscope magnification. The minimum distance that 
can be measured is I pixel. For estimation of the particle 
section sizes with a measurement error of less than 
10 pet, measurements have to be carried out at a micro- 
scope magnification where the particle image sizes are 
of the order of 10 to 20 or more pixels, so that the mea- 
sured sizes have an error of less than 1 pixel (i.e. y less 


than 10 pet). At such magnifications, for a number of 
important classes of microstructures, on the average about 
10 to 50 particles are present in an image frame (field 
of view). Further, in most of the cases, microstructural 
distances up to about 5 to 10 times the mean free path 
(average uninterrupted surface to surface distance be- 
tween the particles through the matrix) have to be mea- 
sured to obtain meaningful distance-distribution data. 
Thus, it is necessary to measure the distances between 
the particles that may not be in the same image frame 
or field of view (Figure 1). To measure such distances, 
it is obviously necessary to know the distances between 
the different image frames on the metallographic plane 
with an error of less than few r pixels. Alternately, a pro- 
cedure must be developed to obtain the ( X , Y) particle 
center coordinates of the particles in all the different im- 
age frames referred to the same (0. 0) origin ; the dis- 
tance between the centers of the particles situated in 
different image frames can be then computed by using 
simple distance formula of coordinate geometry. This 
approach is adopted in the present work. The limitation 
on the size of the image frame (typically 512 x 5 1 2 pix- 
els) also creates another problem due to the associated 
“edge effect. " It is necessary’ to properly account for the 
particles that are only partly in the image frame, and 
their centers may be outside the image frame. Ignoring 
the particles on the edges of the image frame may in- 
troduce a bias of unknown extent in the measured dis- 
tance distribution functions. This bias due to edge effect 
can be eliminated for all the practical purposes if a large 
number of contiguous fields of view (say about 35 to 
50) can be precisely “pasted ” together by developing an 
appropriate image analysis software. In the present work, 
such software is developed for the VIDAS image anal- 
ysis system supplied by Carl Zeiss, Inc. (Germany). 
However, the basic algorithm is applicable to any mod- 
em image analysis system. Once different contiguous 
image frames are pasted together, the distance between 
the consecutive frames is simply equal to length of the 
image frame ( i.e 512 pixels in the present case), and 
hence the (X y Y) coordinates of the particle centers in the 
different image frames can be referred to the same 



Field of view 1 


Fig I — Measurement of distances between particles located in dif- 
ferent fields of view Note that to calculate distance from (A\ K) to 
<A". >"), it is necessary to have the same reference origin for both 
fields of view 
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„ o) origin by simple translation of the coordinates. 

3 The resulting coordinates and the respective panicle sizes 

4 then form the basic raw data necessary to compute any 
'» distance distribution. The algorithm for the image anal- 

* ysis software to paste the image frames and to extract 

* the (X, Y) center coordinates and the particle image size 
It data is as follows. 

The first image frame (field of view) is selected ar- 
bitrarily and stored in the memory of the image analyz- 
es ing computer. The right border (having about a 50-pixel 
^ width) of this image is pasted on the left edge of the 
, next approximately contiguous live-image frame (field 
of view), and the microscope stage is moved and ad- 
i, j us ted so as to fit the right border of the previous frame 
to the left edge of the live image within about 10 to 20 

- pixels; this rough matching is done manually by the op- 
» erator. At this point, the image of the live frame is grabbed 

and translated pixel by pixel until it perfectly matches 
with the right border of the previous field displayed on 
the screen.* Once a satisfactory match (within one pixel) 

•Note that the image can be translated by a preset number of pixels 
via a combination of commands available in the image analyzer 

is achieved, the second image frame is stored The same 
procedure is continued with successive image frames. In 
this manner, a composite global frame of precisely con- 
tiguous 35 image fields (5x7 field of view) is stored 
in the computer; the upper limit on the number of image 
fields is set by the hard disk size and not by the image- 
analysis procedure. The successive image fields in the 

* composite global frame are numbered in a serial order 
The top left-hand comer of the first frame is designated 
as (0, 0) and is origin of the coordinate system (Figure 2). 

Although the composite global image consisting of 35 
fields of view is stored in the computer memory, the 
image analyzer can characterize only one image frame 
of the size 512 x 512 pixels at a time, and hence the 
creation of a composite global image by itself does not 
eliminate the edge effect. To solve this problem, a mea- 

- surement frame of one-fourth of the size of 512 x 512 
pixel frame is defined. The composite global image is 
divided into a number of such frames. For example, in 
Figure 3, the composite image of 16 fields of view is 



Fig. 2 — Coordinate reference for a global image made of 16 contig- 
uous fields of view 


divided into 49 contiguous measurement frames. One* 
fourth of each of these stored image frames, containing 
the corresponding part of the measurement frame, is re- 
called on the analysis screen of the image analyzer; the 
measurement frame is at the center of the screen 
Figure 4 shows this arrangement. In Figure 4, the dashed 
lines are the boundaries of the measurement frames, the 
measurements are performed on the measurement frame 
located in the center and completely inside the analysis 
screen. The measurements are performed only on those 
particles whose centroids are inside the measurement 
frame (for example, dark particles in Figure 4) The image 
analyzer gives (X', Y f ) centroid coordinates and particle 
sizes of these selected particles. At these stage, the cen- 
troid coordinates (X', Y') are with respect to the 
(O', O') origin at the bottom left-hand comer of the 
measurement frame , whose location with respect to 
the fixed origin (0, 0) of the composite global frame 
(Figure 2) is precisely known. Thus, (X, ) ) centroid co- 
ordinates with respect to the fixed (0, 0) origin of the 
composite global frame can be calculated from local 
(X', Y') centroid coordinates by simple coordinate trans- 
lation. The (X, K) centroid coordinates and correspond- 
ing particle sizes are stored in a File as a series of triplets 



Fig. 3 — Division of composite global image frame made up ot 16 
contiguous fields of view (solid lines) into a set of measurement frames 
(dashed lines) 



Fig 4 — The measurement frame for image analysis 
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of numbers. The particles, which are partly in the mea- 
surement frame and whose centroids are outside the 
measurement frame under consideration (white particles 
in Figure 4), are accounted for and measured when the 
corresponding measurement frames containing their re- 
spective centroids are analyzed.* In this manner, cen- 

*These particles would be excluded from measurements if the mea- 
surements were performed on random disconnected fields of view; 
this is the edge effect. Due to this edge effect, the microstructural 
distances involving these particles would not be included in measured 
distance distributions, resulting in a bias of unknown extent 

troid coordinates and particle size of different particles 
are measured, and they are reported only once. There is 
a band of exclusion of about 100-pixel width all around 
the composite global image (consisting of 35 fields of 
view), which does not belong to any measurement frame 
(region bounded by broken lines and global frame 
boundaries in Figure 3), and hence no measurements are 
performed on the particles whose centroids are in this 
band. The error due to this “global” edge effect is ex- 
pected to be very small, and it decreases with the in- 
crease in the number of fields of view in the composite 
global image. 

The data on the ( X , T) centroid coordinates and the 
particle section sizes form the input for the computer 
program used to calculate the K function, radial distri- 
bution function, and the nearest-neighbor distribution 
function pertaining to the 2-D metallographic section on 
which the measurements have been performed. To cal- 
culate the K function, the program proceeds as follows 
Let m be the total number of particles on which the 
measurements are performed, and let (X,, T,), (X 2 , Y?). 

. . . (X,, T,), . . . (X m , Y m ) be the centroid coordinates 
of the particles having the sizes r,, r 2 , . . . r,, . . . r m . 
The program calculates the distances between every par- 
ticle (X,, Y,) and all other (m- 1) particles. It then com- 
putes the number of particle centroids at a distance less 
than r from the particle at (X,, Y ’), and repeats this cal- 
culation for all values of the index /, from 1 to m ; the 
average value of this quantity (properly normalized bv 
magnification) is the estimate of the K function for that 
value of r. The procedure is repeated for all the desired 
values of r to obtain the estimate of the function k(r) for 
the 2-D metallographic section. The program also com- 
putes the total number of particle sections per unit area 
of metallographic plane, jV a . These estimates of k(r) and 
N a together with Eq. [4] are then utilized to compute the 
radial distribution function g(r). The pair correlation 
function can be also estimated by using a similar pro- 
cedure. The nearest-neighbor distribution can be esti- 
mated simply by recording the distances between the 
particle centroids and respective nearest panicles. 

IV. A PRACTICAL EXAMPLE 

The electrical conductivity and other physical prop- 
erties of polymer matrix composites containing hollow- 
spherical carbon particles in a polymeric resin matrix are 
of interest for certain electrical shielding applications. In 
such a material, the conductivity of the polymeric matrix 
is significantly lower than that of the carbon particulate 
Thus, the conductivity and other electrical properties of 


the composite are very sensitive to the spatial arrange- 
ment of the carbon particles in the microstructure and 
the interparticle contacts among these particles. The ex- 
perimental measurements of the microstructural distance 
distributions of the hollow carbon particles in the poly- 
mer resin matrix have been performed to demonstrate the 
image analysis procedure developed in Section III. The 
measurements are performed on this composite because 
it is easy to fabricate in the laboratory, its electrical 
properties are sensitive to the microstructural spatial ar- 
rangement, and the material is of practical interest for 
certain aerospace applications. 

A. Materia I 

The hollow spherical carbon particles were supplied 
by Carbosphere Inc. (United Stales). Figure 5 shows the 
size distribution of these particles. The polymeric resin 
(araJdite) and the hardener were supplied by CIBA GEYG1 
(France). The resin, the hardener, and the particulate 
powder were mixed in the appropriate proportions, and 
the mixture was heated in vacuum to degas the material 
and to decrease the viscosity of the polymer. It takes 
about 10 hours for the polymer to harden completely, 
and hence to minimize the particle segregation due to 
gravitational effect (the specific gravity of particles is 
significantly lower than the liquid polymer), the molds 
were slowly rotated during the hardening of the polymer. 
The specimens having particulate volume fraction of 0.29 
and 0.47 were prepared in this manner. 

B. Materiallography and image Analysis 

The specimens were cut, mounted, and polished on a 
series of SiC polishing papers, followed by polishing on 
diamond cloths using automatic polishing wheels. 
Figure 6 show s a typical microstructure of the composite 
observed on a random metallographic plane. The carbon 
particles are round in shape, and there is a range of par- 
ticle sizes. In Figure 6, the dark regions are the “hol- 
lows” created by the intersections of hollow spherical 
particles with the sectioning plane. The bright regions 



Sire (jjjh) ► 

Fig. 5 — Frequency distribution of the panicles The volume fre- 
quency is the fraction of total particle volume in a given size range 
The curves given by image analysis and granulometry are in reason- 
able agreement with each other 


metallurgical and materials transactions a 


VOLUME 2bA, JUNE IW - 145? 




Fit 6 Micro’.irui luri's observed in a metalloeraphic section (ui ->d nr, , . , 

ih ' ^ P 1 ' pan kul.no volume Irattion specimen a. magnification 5(Xh,™? r V' ” nk ' ,' raC "' >n specimen at magnification 200 times, 
-m mill's and „/t 47 pa paniculate volume fraction specimen at magnification SlNlir' 1 " " " H friK ’"° n speci,ncn *' m »g"*fic«ion 


are due to near-tangential intersections between the car- 
bon panicle shells and the metallographic plane. The 
contrast between the panicles and the matrix is not good 
and hence a number of image transformations are nec- 
essary to develop the appropriate contrast via image 
processing /-' b 

The image analysis and measurements were per- 
formed at a microscope magnification of S00 times 
yielding a pixel si/e of 0.325 M m. At this magnification.' 
i on the image analy/cr screen, the average panicle size 
’j ,s aboul r ’ t) l 1,xcls - permitting estimation of panicle sizes 
and interpanicle distances with an error of less than 
5 pet In the. low volume fraction specimen (29 pet par- 
ticulate ». the measurements were perfonned on 105 fields 
of new tat magnification 500 times) spread over three 
independent composite global image frames (35 fields of 
view each) In the specimen containing 47 pet carbon 
paniculate, the measurements were performed on 70 fields 
ol view at magnification 500 times spread over two in- 
dependent composite global image frames, each com- 
posed ol 35 lields of view . The image analvsis procedure 
was used to obtain the iX. V) centroid coordinates and 
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the panicle section sizes, as described in the last section 
The kin function, the radial distribution function g(r ), 
and the nearest-neighbor distribution function p{r) were 
estimated Irom these data by using a computer program 
written m the C language. These results are discussed 
subsequently 

C. Results iitu I Discussion 

Figure 7 shows the variation of the function k(r) with 
the nucrostmctural distance r between the panicle sec- 
tion centroids observed in a metallographic plane. Recall 
(hat kit t represents an average number of panicle section 
centroids at a distance less than r from an arbitrary' par- 
ticle centroid The value of k(r) increases monotonicaily 
with /. n approaches 0 as r approaches 0, and n tends 
to infinity as / goes to infinity, as expected. At any given 
value of / . £<r) is higher lor the high volume fraction 
specimen than lor the lower volume fraction specimen. 
Figure H shows the radial distribution function g(r) for 
the two spes miens plotted against the normalized micro- 
structural distance (r/r). where f is the average particle 
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Fig 7 — Estimated K-function of the carbon panicle sections ob- 
served in a metaHographic plane. 



r/t y 

Fig. 8 — Estimated radial distribution function of the carbon particle 
sections observed in a metaHographic plane 


section size. The first maximum in the radial distribution 
function occurs at about (r/f) equal to 3.0 for the high 
volume fraction specimen and at (r/f) equal to 3.5 for 
the low volume fraction specimen. These distances in- 
dicate the maximum likelihood of finding a particle sec- 
tion centroid from any given particle in a metaHographic 
plane. It is important to recognize that the normalized 
radial distribution functions are different for the two 



Fig. 9 — Estimated nearest-neighbor probability distribution of the 
distances between carbon particles in the polymer matrix composite. 
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specimens having different volume fractions. This is in- 
dicative of the intrinsic effect of volume fraction on the 
spatial arrangement of the particulate. The radial distri- 
bution function approaches 1 at distances larger than about 
five times the average particle section size (although very 
small perturbations around this mean value exist). Recall 
that for randomly distributed particle centroids, g(r) is 
equal to 1. This shows that for the present specimens, 
particle-particle affinities/repulsions are negligible at 
distances larger than about five times the average par- 
ticle size. This information is useful for the computer 
simulation models. Due to memory limitations, in the 
computer simulation studies, only a finite extent of 
microstructural volume can be simulated, and it is of in- 
terest to simulate the minimum microstructural volume 
(called representative volume) that gives the reliable 
process information. For the processes where only short- 
range particle-particle interactions are important, the dis- 
tance beyond which g(r) approaches 1 should provide a 
realistic measure of the “representative volume element” 
for computer simulation models. 

Figure 9 shows the estimated nearest-neighbor distri- 
bution function for the two specimens under investiga- 
tion, plotted against the normalized nearest-neighbor 
distance. The average nearest-neighbor distances are 17.0 
and 14.6 fxm for the specimens having 0.29 and 0.47 
particulate volume fraction, respectively . 
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